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We study the Josephson effect in a triplet superconductor-ferromagnet-singlet superconductor 
junction. We show that the interaction of tunneling Cooper pairs with the interface magnetization 
can permit a Josephson current at the lowest order of a tunneling Hamiltonian perturbation theory. 
Two conditions must be satisfied for this to occur: the magnetization of the ferromagnet has a 
component parallel to the d-vector of the triplet superconductor, and the gaps of the superconductors 
have the same parity with respect to the interface momentum. The resulting charge current displays 
an unconventional dependence on the orientation of the magnetic moment and the phase difference. 
This is accompanied by a phase-dependent spin current in the triplet superconductor, while a 
phase-independent spin current is always present. The tunneling perturbation theory predictions are 
confirmed using a numerical Green's function method. An analytical treatment of a one-dimensional 
junction demonstrates that our conclusions are robust far away from the tunneling regime, and 
reveals signatures of the unconventional Josephson effect in the critical currents. 

PACS numbers: 74.50.+r, 74.20.Rp 



I. INTRODUCTION 

The physics of spin-triplet superconductors (TSs) is 
much richer than their spin-singlet superconductor (SS) 
counterparts due to the spin degree of freedom of a triplet 
Cooper pair. Although a bewildering variety of triplet 
pairing states are in principle possible for any crystal 
symmetry^ there are only a handful of systems where 
a TS state has been well established, e.g. UPt3 and 
S^RuC^P^Even in these cases, many questions re main 
about the precise form of the TS order parameterP^ 
Much effort has therefore been directed at developing 
experimental tests capable of unambiguously identifying 
TS. A promising route is to incorporate the candidate 
TS into a heterostructure device and search for signa- 
tures of t he o dd-parity orbital state in tunneling mea- 
surements!^ Alternatively, the spin part of the Cooper 
pair wavefunction can be probed by bringing the TS into 
contact with a ferromagnet (FM). This is manifested by 
the crucial role of the relative orientation between the 
vector order parameters of the TS and the FM, the d- 
vector and magnetization respectively, in controlling the 
physics of the deviceP^ 

Theoretical investigations of the Josephson effect be- 
tween TS have revealed many remarkable consequences 
of their intrinsic spin structure. For example, the 
spin of the Cooper pair permits the existence of a 
Josephson spin current between TS with misaligned d- 
vectors p 4 * 15 * simil ar to the spin supercurrent observed 
in superfluid 3 He! 16 * 17 ^ Another notable proposal is the 
triplet superconductor-ferromagnet-triplet superconduc- 
tor (TFT) junctionp^HIH where the coupling of the 
Cooper pairs' spin with the exchange field in the barrier 
causes a sign reversal of the current as the orientation of 



the exchange field with respect to the d-vectors is var- 
ied. This is in stark contrast to the well-known 0-7T tran- 
sition in magnetic junctions between SSsp2H24] which is 
independent of the orientation of the barrier magnetiza- 
tion. Such an anomalous 0-7T transition would therefore 
be extremely strong evidence of the triplet state of the 
superconductors. On the other hand, creating a TFT 
junction from any of the TS candidate materials is ex- 
perimentally challenging due to the high purity require- 
ments needed for the TS. Even with recent success in 
growing superconducting thin films of S^RuCUpSl such 
devices will likely remain hypothetical for some time, due 
to the greater challenges posed by growing a layered het- 
erostructure. More immediately plausible is to create a 
triplet superconductor-ferromagnet-singlet superconduc- 
tor (TFS) junction by coating the S^RuC^ thin film with 
a ferromagnetic layer and then contacting to a conven- 
tional superconductor. 

For junctions between superconductors of like parity, 
the lowest harmonic in the Josephson current vs phase 
difference <f> relation is usually sin(0), which originates 
from tunneling processes involving only a single Cooper 
pair. This term is necessarily absent in a nonmagnetic 
Josephson junction between an SS and a TS, however, 
due to the opposite parity of the pairing statespS Instead, 
the parity conversion can only occur in processes involv- 
ing the coherent tunneling of even numbers of Cooper 
pairs, and so sin(2(/>) is the leading harmonic in the cur- 
rent vs phase relationP^^The coupling between a single 
tunneling Cooper pair and the magnetic degrees of free- 
dom in a magnetically-active barrier, on the other hand, 
can accomplish the conversion between singlet and triplet 
spin states,^ hence generating a lowest-order Josephson 
coupling in the sense discussed above. An example of 
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such a magnetic interaction is the intrinsic spin-orbit cou- 
pling expected to occur at the junction interfaceP^H 
This has been proposed as the origin of the unexpectedly 
large Josephson currents injunctions between single crys- 
tals of Sr2Ru04 and conventional s-wave SS, and the pro- 
nounced dependence on the crystal face upon which the 
Josephson contact is madepSHSU i n contrast, relatively 
little work has been done for a ferromagnetic tunneling 
barrier. Previous studies have found a highly unconven- 
tional Josephson charge current in a TFS junction which 
is even in the phase difference and odd in the compo- 
nent of the magnetization parallel to the d- vector of the 
TSpSHUl The origin of this current, and the conditions 
under which it occurs, nevertheless remain obscure. 



triplet ferromagnetic singlet 

superconductor barrier superconductor 




FIG. 1. (color online) Schematic diagram of the Josephson 
junction considered in this paper. The d-vector of the TS de- 
fines the x-axis, and we restrict the magnetization to the x-y 
plane, making the angle a to the x axis. The phase difference 
between the TS and SS is given by <j>. 



It is the purpose of this paper to perform a detailed 
theoretical study of the Josephson effect in a TFS junc- 
tion, shown schematically in Fig. ([!]). Similar to the 
TFT junction, we find highly unconventional Josephson 
physics which originates from the coupling of the barrier 
moment to the spin of the tunneling quasiparticles com- 
prising the Cooper pairs. Using a tunneling Hamiltonian 
perturbation theory, we show that there is a lowest-order 
Josephson charge current when the orbital pairing states 
of the superconductors have the same parity with respect 
to the interface momentum, and the magnetization has 
a component parallel to the d-vector of the TS. There is 
also a spin current in the TS, which has both a phase- 
dependent and a phase-independent contribution. The 
latter is a universal spin supercurrent which appears at 
triplet superconductor-ferromagnet interfaces and is due 
to spin-dependent reflection processes.^ We test the pre- 
dictions of the perturbation theory using both lattice and 
continuum models of the junction. In the lattice the- 
ory we survey a wide selection of different orbital pairing 
states in the SS and TS. The focus of the continuum 
theory, on the other hand, is to understand the role of 
resonant tunneling through the Andreev bound states at 
the junction interface. Both approaches yield excellent 
agreement with the perturbative analysis at sufficiently 
high temperatures. Although deviations from pertur- 
bation theory become more severe as the temperature 
is lowered, it nevertheless remains qualitatively correct 
down to zero temperature. 



This paper is organized as follows. In Sec.[TT]we present 
a perturbative theory for the Josephson effect in the 
TFS. The predictions of this section are confirmed first 
in Sec. |III| by numerically determining the currents in a 
microscopic lattice model of the junction, and then by 
analytical calculation for a continuum model in Sec. |IV| 
Our concluding discussion is given in Sec. [V] 



II. PERTURBATION THEORY 
A. Hamiltonian of the junction 

The Hamiltonian of the TFS junction is written 

H=H T s+Hss+H t un+H Icf . (1) 

Here Hts and Hss describe the bulk triplet and singlet 
superconductors on either side of the barrier. We have 

k / 

(3) 

where 9 vM = (<v )kff ,c l/4a ,4,_ k)t ,4 1 -k P |.) T and 'W 
(4 k ) are the fermion annihilation (creation) operators 
for a spin-cr quasiparticle with momentum k in the v = s 
(t) singlet (triplet) superconductor. The bare dispersion 
in each superconductor is e^k- We only consider unitary 
equal-spin-pairing states for the TS. Without loss of gen- 
erality, it is convenient to let the vector order parameter 
of the TS define the x-axis, i.e. d k = At t ke x . Other 
orientations of d k can be achieved by spin rotation of 
the system, and do not result in new physics. The global 
phase difference between the TS and the SS is <j>. 

The TS and SS are connected by the tunneling Hamil- 
tonian, which we define 

%un= 13Z] T ^k'4,k, CT c W ,k',a' , ( 4 ) 

u=s,t k,k' a,a' 

where V = s(t) for v = t(s). In order to properly 
model the interaction of the tunneling quasiparticles with 
the magnetic moment of the barrier we require explic- 
itly spin-dependent tunneling matrix elements. Follow- 
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ing Ref. [201 we also include a reflection Hamiltonian 

Hrcf = Z Z Z ^,k,k'4,k )CT c ^k' -a > (5) 

to properly account for the interaction of the quasipar- 
ticlcs with the FM layer. We only include spin-flip re- 
flection processes in "H re f, since spin-preserving reflection 
processes clearly do not contribute to either the charge 
or spin currents. 

1. Ansatz for tunneling and reflection matrix elements 

Our perturbation analysis crucially relies upon the 
form of the tunneling and reflection matrix elements. In 
particular, it is necessary to include the phase shift ac- 
quired by the quasiparticles during the various tunneling 
or reflection processes.^ Although it is in principle pos- 
sible to determine the tunneling and reflection matrix 
elements from a more fundamental Hamiltonian, here we 
motivate a phenomenological form by comparison with 
an exactly-solvable scattering problem. 

A common model for the tun neling barrier is a 8- 
function potential™ ™™! i.e V(r) = U a S[z) + 
Um • o5 (z), where Uq is the charge scattering potential 
and XJm = Um (cosae^ + sinae y ) is the magnetic scat- 
tering potential. It is straightforward to evaluate the 
scattering matrix for this potential, and hence determine 
the transmission and reflection coefficients, t v ^, a i (k, k') 
and r vaa i (k, k') respectively. In particular, the trans- 
mission and spin-flip reflection coefficients are both small 
when the magnetic potential is large, i.e. the dimen- 
sionless parameter g = kp\U mI/^Ep 3> 1, where kp is 
the Fermi wavevector and Ep is the Fermi energy. It is 
reasonable to expect that in this limit we should have 

T „£.y ~ W'(k,k') and ~ r„, CT ,_ CT (k,k'). By 

requiring that a tunneling or reflected quasiparticle ac- 
quires the same phase as in the exact solution, we hence 
have the following ansatz for the matrix elements 

T%,v = ^T 8p (kz,K)K*\r (6a) 

iup iaa 

T U£y = —T sf {k z X)KK ' (6b) 

iup iaa 

R^^—R sf (k z ,K)S khk>iil (6c) 

where v = — 1 (+1) as a factor for v — t (s) as a subscript. 
Crucially for our analysis, the phase shift acquired dur- 
ing spin-flip tunneling or reflection depends on the initial 
spin a and the angle a of the magnetic moment in the 
x-y plane [see Fig. Q]. We assume that the T sp (k z ,k' z ), 
T s f {k z ,k' z ), and R s ^[k z ,k z ) are real functions, indepen- 
dent of g, and satisfy 

T s P(k z ,k' z ) = T sp (-k Zl -k' z ) = T s P(k' z , k z ), (7a) 
T s f(k Z7 k' z ) = -T s f(-k z , -k' z ) = T s f(k' z , k z ) , (7b) 



R sf (k*,K) = -R sf {-k z ,-k' z ) = -R s f{k' z ,k z ). (7c) 

These conditions originate from both the comparison to 
the scattering coefficients and the requirement that %tun 
and "Href be Hermitian. Note that the (5k y in Eq. |6| en- 
sures the conservation of momentum parallel to the bar- 
rier, a consequence of the translational invariance along 
the interface™ 

Although we have motivated our ansatz Eq. (|6| by 
comparison with the 5-function barrier, we expect our 
approach to be of more general validitypD In particular, 
the spin-dependent phase shifts acquired by the tunneling 
quasiparticles should be robust to other choices of barrier 
model, since they are a consequence of the orientation of 
the magnetization. 



B. Perturbation theory 

The number operator for particles in each spin sec- 
tor of the two superconductors is given by N Vt(J = 

c l k o c v k a- F rom this we define the associated par- 
ticle currents 

r a = -v{d t N v ^) . (8) 

We proceed by expanding the S matrix to lowest order 
in %' = Htun + "Href, hence treating the tunneling and 
reflection processe s as a perturbation of the Hamiltonian 
"Ho = n T s + "HssP 21421 This is justified so long as the 
tunneling and reflection matrix elements are small, which 
by our ansatz Eq. ^ holds in the g 1 limit. By the 
Kubo formula, we then have 

F a = -iv f dt'([d t N^(t),H'(t')}), (9) 

where the time-dependence is given within the interac- 
tion picture, i.e. 0(t) = e mot Oe- mot , with H = 
Kq + l^v J2a N v ,<t and Kq is the associated grand 
canonical Hamiltonian. Since we are only interested in 
the DC Josephson effect, we take the same chemical po- 
tential in both the TS and SS, i.e. fj, s = fj>t = /i- We 
hence find for the terms in the commutator in Eq. Q 

d t N Via (t) = i {B-^(t) - BZ'-"{t)} 

+i J2{Ar(t)-A^(t)}, (10) 

= Z Z w (*) + Z Z & (*) . ( n ) 

where we introduce the operators 

Al*\t) = ^T^4^(t)c vM , t ,,(t), (12) 

k,k' 

B^\t) = £ Rli >k 4^(t)c uXt ,,(t) . (13) 

k,k' 
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The time-dependence of the fermion operators in these 
expressions is given by c v k ff (i) = e lKot c v k a e~ lKot . 

Following standard arguments j 15 * 20 * 4 ^ we write the cur- 
rent Eq. ([9} as 

II = 2vlm { (oj = 0) + (u = 0) } , (14) 

where the retarded correlation functions and 
ty T J% (ui) give the contributions from tunneling and reflec- 
tion processes, respectively. They are obtained by ana- 
lytic continuation ioj n — > cu + i0 + of the corresponding 
Matsubara functions 

^,,K)= f dre^ T {T t A^(t)A^"{Q)){^) 
Jo 



*uA">n)= dre w - T (T T B-^(r)B-^(0)). (16) 
Jo 

The Matsubara functions are evaluated by using Wick's 
theorem to expand the two-particle correlators. We 
hence find the particle currents in the TS 

4 = -^£[tf s/ («)]X*t, 

k,k' 

xl m ^-^l^(k,k') 

-C't.k-C't.k' 



J2T°P(k z ,k' z )T°f(k z ,k' z )5^, 



k,k' 



;B e{e«(^ y k ^ k, U. |t(kik >), (17) 



and the SS 



II = 



2 cos (a) 
xRe<! e ! 



k,k' 

A:kA/k ' ' K,<k.k': 



^s.k^t.k 

Here we utilize the convenient short-hand notation^ 
F vvl {^)= f{EvM) ~ f{Eu ' M ' ] 



(18) 



where 



E vM = VKk-^) 2 + |A,, k | 2 



(19) 



(20) 



is the dispersion in the v superconductor and f(E) is 
the Fermi distribution function. Equations (JT7J) and ( fl8| ) 
are the central results of our perturbation theory anal- 
ysis, as they give all lowest-order contributions to the 
particle current in each superconductor. Note that in 
the TS we have a contribution from reflection processes, 
whereas in the SS only tunneling processes contribute to 
the current. The lowest-order (i.e. single-Cooper-pair) 



tunneling processes involve both a spin-preserving and a 
spin- flip tunneling event, which is necessary to transform 
the spin singlet Cooper pairs of the SS into the S z = ±h 
triplet Cooper pairs of the TS and vice versa. The par- 
ticle current in the TS [Eq. (17)] depends on the spin 
a through the spin-dependent phase shifts acquired dur- 
ing the scattering; in contrast, the particle current in the 
SS [Eq. ( 18 )] is the same for each spin orientation, as 
required by the singlet pairing state. 

The important role of spin-flip tunneling implies a 
transfer of spin to the FM tunneling barrier. In our cal- 
culation, however, we regard the magnetic moment of 
the FM barrier to have fixed magnitude and direction. 
Including the response of the FM to the injected spin 
current is a challenging problem, requiring a noncquilib- 
rium treatment that is beyond the scope of the current 
manuscript. 



C. Charge and spin currents 



We can use either Eq. ( 17 1 or Eq. ( 18 ) to calculate the 
Josephson charge current I c = — e{K + IY), as by charge 
conservation this is the same in each superconductor. We 
hence find 



I, 



k.k' 



^^#4^K(k,k'). 

I Ai.k-ftt.k' J 



(21) 



The charge current strongly depends upon the orienta- 
tion of the magnetic moment through the cos(a) factor. 
This implies that reversing the direction of the barrier 
moment also reverses the sign of the current, as was pre- 
viously observed in Ref. HH1 The origin of this factor 
is the interference of the particle currents in each spin 
sector of the TS, which are phase-shifted with respect 
to one-another by ±2a as a consequence of the spin-flip 
tunneling, see the second term in Eq. (17 1. 



We extract the current vs phase relationship by ex- 
amining the summand in Eq. (21). In order to have a 



lowest-order Josephson effect, we require that the prod- 
uct T s f(fc z ,fc;)T^(^,^)A; k A t ^ ( 5 k|| ^ | not be odd in 

any component of k or k'. From Eq. ^ we hence deduce 
that the gaps A* k and A tk / have opposite parity with 
respect to the z-componcnt of the wavevector; equiva- 
lently, the gaps must have the same parity with respect 
to the interface momentum. For an s-wave SS, therefore, 
there is a Josephson current 



I c cx — cos(a) cos(0) , 



(22) 



when the TS has p z -w&ve symmetry. The proportional- 
ity constant is determined by the details of the junction, 
such as the normal-state dispersion or the structure of the 



interface. On the other hand, Eq. (21) is vanishing for 
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TS gap symmetry 
Pz, Pz +ip v 
Py 
Pz 
Py 

Pz + tPy 



SS gap symmetry 

s, d„2 __2 



s, d„ 



charge current [Eq. (21 1 
I c cos(a) cos(0) 



I c cos(a) cos(0) 
I c cos(a) sm(4>) 



spin current [Eq. ( 23 1] 



I Si r sin(2a) + sin(a) sin(0) 

I s , r sin(2a) 

I s , r sin(2a) 
I S: r sin(2a) + I a j sin(a) sin(</>) 
I S) r sin(2a) + I S: t sin(a) cos(0) 



TABLE I. Table showing the Josephson charge and spin currents appearing in the lowest order of our perturbation theory for 
different combinations of TS and SS gap symmetries. In the interests of brevity we restrict ourselves to gaps lying in the y-z 
plane. The Josephson current amplitudes are expected to satisfy I c ~ g~ 3 , I s ,r ~ g~ 2 , I s ,t ~ g~ 3 in the tunneling limit g 2> 1. 
The numerical values of these terms depend upon the details of the junction, e.g. the normal-state dispersion, the structure 
factor of the gaps, the properties of the tunneling region, etc. 



a Py-wave TS, and so a Josephson current only appears 
in the next order of perturbation theory. The current vs 
phase relation is then I c oc sin(2<^>), as for the nonmag- 
netic barrier™ A list of the Josephson charge current 
vs phase relationships predicted by Eq. (21) for different 
combinations of orbital pairing states is given in Table [I] 
The spin current is only present in the TS, as the sin- 
glet Cooper pairs do not carry spin. The spin current is 
polarized along the z-axis and is given by 

X (4 -n) 

k.k' 



I, 



xRe 



2 sin(a) 



F,,fk,k') 



9° 



E 

k,k' 



T s P(k Z7 k' z )T^(k z ,k' z )S: 



kn,k' 



xlm 



,^ A s,k A «,k' 



^,t(k,k') 



(23) 



The first term in Eq. ( 23 ) is the spin current due to spin- 



flip reflection from the interface, where a spin S z = ±H 
Cooper pair is reflected as a S z = Cooper pair! 11 | 2 ° ^ 
As this involves two spin-flip reflection events, the Cooper 
pair acquires a phase shift of ±2a and the magnitude of 
this term goes as ~ g~ 2 . Furthermore, due to its ori- 
gin in reflection processes, it is sensitive to the orbital 
structure of the triplet gap, e.g. the sign reversal of the 
p 2 -wave gap upon specular reflection gives the Cooper 
pairs an additional 7r phase shift relative to the p y -wave 
case, and there is hence a sign difference in the respec- 
tive reflection spin currents. In contrast, the second term 



in Eq. ( 23 1 originates from the interference of the spin-t 



and -J. tunneling particle currents, similar to the charge 
current, and goes as ~ g~~ 3 . We therefore expect that the 
reflection spin current dominates the tunneling spin cur- 
rent in the tunneling limit g 3> 1. Whereas the reflection 
spin current is always present, the condition for the tun- 
neling term in Eq. ( 23 1 to be non-zero is the same as for 
a finite second-order charge current. Again considering 
the case of an s-wave SS, we thus find that for a p z -w&~ve 



n=1 2 



L L+1 




TS 



SS 




FIG. 2. (color online) Schematic diagram of the two- 
dimensional lattice model of the TFS junction. The length 
of the ferromagnetic layer and the junction width are L and 
M unit cells, respectively. The on-site magnetic potential 
UmCO m t ne ferromagnetic layer is oriented in the x-y plane. 



TS there is the spin current 

7 1 
Is z oc sin(2a) + sin(a) sin(^) , (24) 

9 9 

where 7 is a numerical constant; for a p y -wave TS, in 
contrast, only the first term is present, and has opposite 
sign. The Josephson spin currents are given in Table [T] 
for various combinations of orbital pairing states. 



III. LATTICE MODEL OF THE JUNCTION 

In this section we test the predictions of the perturba- 
tion theory by calculating the charge and spin Josephson 
currents in the tunneling regime of a two-dimensional mi- 
croscopic lattice model of the TFS junction. Using the 
recursive Green's function methodp^ we examine repre- 
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sentative examples of the different orbital combinations 
listed in Table Q] 



A. Formulation 

We write the Hamiltonian describing the TFS junction 
in real space as 



(25) 

where Vl'(r) = (ijj f (r), (r), V4( r )' (r))' r is the vector 



of field operators and the sum in Eq. (251 is over the 



lattice sites. For the normal state Hamiltonian Ho(r,r') 
we assume a two-dimensional square lattice tight-binding 
model in the y-z plane 



£o(r,r') 



-t{6 r 
-jj,5 TT 



00 



U r,r — z 

U M (r) 



&S r r 



-y} 
(26) 



The vectors are represented as r = nz + my where z and 
y are the unit vectors of the tight-binding lattice in the 
z and y directions, respectively. In the y-direction, we 
apply the periodic boundary condition. The exchange 
potential Ujvf(r) is only non-zero in the ferromagnetic 
barrier region of the junction. A sketch of the lattice 
model is shown in Fig. |2]). The pair potentials for the 
different pairing symmetries considered here have the fol- 
lowing real-space forms: 



A(r,r') = < 



Ae l ^8 r yia y 
(Ae-^/2)ia y [S r y +z+y 

i(A/2)a z [£ r ,r'+z - <5r,r'-z 

i(A/2)a z [S r y +y - c> r:r '_ y 

(A/2)<7 2 [iS r y +z - iS r y-, 



^r,r' — z— y 



3r,r'+y 



J r,r'+z— y 



<5r,r'-yj 



<5r,r'— z+y] 



s-wave, 

d Z y-w&ve, 

p z -wave, 

Pj,-wave, 

Pz + ipy-wave 



(27) 



The pairing potentials for the TS assume that d is di- 
rected along the x axis in spin space. 

The charge and the spin densities in the n-ih column 
along the z direction are defined by 



M 

p{n) =e E i>U r )i><*( T ) . 

m— 1 a 

h M 

s ( n ) = 2 E E ^( r )*«,^/5( r ) • 

m=l a. ft 



From the equations of motion, 



9 t p(n) =-[H,p{n)], 
d t s(n) =-[H,s(n)], 
we derive the current conservation laws 



(28) 
(29) 

(30) 
(31) 



d t p(n) = -I ce {n) + I ce (n - 1) - S cd (n) , (32) 
d t s(n) = -I se (n) + I se (n - 1) - S sd (n) - S„(n) .(33) 

On the right hand side of these equations we have terms 
that can be interpreted as the divergence of a current. 
The first of these terms are the so-called "kinetic" cur- 
rents which originate from the commutator of the densi- 
ties with the hopping Hamiltonian: 



-^(r)Va(r + z)] 



Ise(n) 



it 



M 



2 EE [^«( r + z )° a ,cs^[3(r) 

m—1 a, 8 



-^(r)ff a ,^/s(r + z)] 



(34) 



(35) 



The remaining terms in Eq. ( 32 1 and Eq. (|33j) are the 
source terms. Specifically S c d[n) is the source term for 
the electric current due to the pair potential, S s d(n) is the 
source term for the spin current due to the pair potential, 
and S.„(n) is the source term for the spin current due to 
the exchange potential. They are represented by 



S c d{n) = 



IT 



EEE[^( r )^( r ')A^(r,r') 

m—1 a, (5 r' 

+V^(r)^(r')A; !Q (r,r')l , 



(36) 



M 



S sd (n) = y E E E [^(r)V>i(r')AA,Mr,r'K, Q 

m=l a,/3,X r' 

+Va(r)^A(r')Al /3 (r,r')<r^ a ] , (37) 

M 

S„(n) = - E E [^W i U MW >< °} a ,pM*)] (38) 

m—1 a, ft 

Currents due to the source terms in spin-triplet super- 
conductor are defined 



I c d{n) = - Y S cd(i) 

n+l<i<no 



(39) 
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hd(n) = - ^2 S ds (i) , 

n+l<i<rio 

I v (n) = - ^2 S„(i), 



(40) 
(41) 



where 



n+l<i<L 



I c (n) = I ce (n) + I cd {n) , 
F° tal (n)=I s (n)+l v (n), 
I s (n) = I se {n) + l sd (n) . 



(44) 
(45) 
(46) 



where no should be in the ferromagnetic layer. We 
are therefore able to re-write the continuity equa- 



tions Eq. (32) and Eq. (33) as 



d t p(n) = -I c (n) + I c (n-l), (42) 
d t s(n) = -I* otal (n) + If tal (n - 1) , (43) 



The averages of these currents are expressed in terms 
of the Matsubara Green's function G defined by 

G(r,r',T-r') = -(T r v]/(r)*t(r')) 

= TY J G{T,T',u n )e- i ^-r') 5 (47) 

where uj n = (2n+ l)irT are the Matsubara frequencies at 
temperature T. Specifically, we write^ 
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G(r,r'. 


l S q(n) 
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(G(r + z 




. M 

= |E T EE Tr 

ni—l uj n r 1 


G(r,r', 


S v (n) 


i M 

= ^ E T E Tr 

m— 1 u;^ 


G(r, r,u n 



A(r',r) 
A*(r',r) 



,[ A(r',r 
^A*(r',r) 

U M (r) x tr 







-U M (r) x <j 




(48) 
(49) 
(50) 
(51) 
(52) 



The recursive Green's function method enables us to nu- 
merically calculate the G reen's function, and hence eval- 
uate the above equations! 44 * 45 ! 



The charge current I c (n) in Eq. (44) is independent of 
n, as required by the charge conservation law. Spin must 
also be conserved, and so we find that I' otal (n) is vanish- 
ing for all n because the spin current cannot flow in the 
SS. This result contradicts the prediction of Sec. [Til that 
there is a Josephson spin current in the TS. The paradox 
can be resolved by noting that neither the perturbation 
theory nor the Green's function method can properly ac- 
count for the transfer of spin to the ferromagnetic bar- 
rier, as in both theories the magnetic moment is assumed 
fixed by the constant exchange potential Ujvf(r). Rig- 
orously accounting for the conservation of spin in such 
a situation naturally leads to the conclusion of vanish- 
ing spin current. It is nevertheless reasonable to identify 
the Josephson spin current with the current l s (n), which 
produces a torque on the ferromagnetic barrier, and to 
hence regard l v (n) as a compensating current necessary 
to maintain the constant exchange potential. Although 
I s (n) depends on n in the ferromagnet, it is independent 
of n in the TS. In the following we only consider this 



spin current in the TS, in order to make contact with the 
perturbation theory. 

In the following we present results for a junction of 
width M = 10 and ferromagnetic barrier length L = 10. 
In units of the transfer integral t we take |Um| = 0.1 
for the exchange potential, and fi — 2 for the chemi- 
cal potential. The pairing potential has weak-coupling 
temperature-dependence, with zero-temperature magni- 
tude Ao = 0.01. We have confirmed that the transport 
properties are qualitatively insensitive to choices of these 
parameters. Since Um x d || z as shown in Fig.[2j the x 
and y components of the spin current are zero. To com- 
pare with the analytical predictions of Sec. [TTJ, we fix the 
temperature at T = 0.5T C in the tunneling regime. 



B. 



-wave singlet superconductor 



In this subsection we present results for the Joseph- 
son currents in an TFS junction between an s-wave sin- 
glet superconductor and each of the three different triplet 
superconducting states listed in Eq. (27). Commencing 
with the p z -wave TS (the p z -F-s junction), in Fig. [3] we 
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FIG. 3. (color online) Josephson currents in the p z -F-s junc- 
tion, (a) Charge and (b) ^-component spin currents as a 
function of <f> f° r fixed a. (c) Charge and (b) ^-component 
spin currents as a function of a for fixed <f>. We choose the 
parameters as M — L — 10, /i = 2t, \TJm\ = O.lt, Ao = O.Oli, 
and T = 0.5T C . 



FIG. 5. (color online) Josephson currents in the (p z + ip y )-F- 
s junction, (a) Charge and (b) z-component spin currents as 
a function of <j> for fixed a. (c) Charge and (b) ^-component 
spin currents as a function of a for fixed <j>. The parameter 
values are fixed as in Fig. [3] 
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FIG. 4. (color online) Josephson currents in the p y -F-s junc- 
tion, (a) Charge and (b) 2-component spin currents as a 
function of <j> for fixed a. (c) Charge and (b) ^-component 
spin currents as a function of a for fixed <j>. The parameter 
values are fixed as in Fig. [3] 



plot the charge and spin currents as functions of the phase 
(f> and angle a. As can be seen in panels (a) and (c), the 
dominant term in the charge current is 



I c = I c cos(</>) cos(a) , 



(53) 



with a much weaker contribution oc sin(20) visible at 
a = 0.57T in panel (a). This is clearly consistent with the 
perturbation theory predictions. The spin current also 
agrees with the perturbative analysis, with the numerical 
results well described by 



I s sin(2a) + I' s sm(<j)) sin(a) 



(54) 



where I s ^> I' s . Indeed, in panel (d) the spin current vs a 
curves at different (f> almost overlap due to the very weak 
^-dependence. The much smaller coefficient of the in- 
dependent term was anticipated in our tunneling Hamil- 
tonian analysis. 

We now consider the results for the p y -wave TS, which 
are shown in Fig. [4] In contrast to the p z -wave TS, the 
dominant contribution to the charge current is oc sm(2<^), 
there is only very weak dependence upon a, and the 
maximum critical current is much smaller. Since the 
oc sin(20) term originates from coherent tunneling of 
two Cooper pairsj^ these results are consistent with 
our prediction of vanishing charge current due to single- 
Cooper-pair tunneling processes. For the spin current we 
find I StZ = I s sin(2a) to excellent approximation, in per- 
fect agreement with the perturbation theory predictions. 
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FIG. 6. (color online) Josephson currents in the p z -F-d yz 
junction, (a) Charge and (b) ^-component spin currents as 
a function of <j> for fixed a. (c) Charge and (b) ^-component 
spin currents as a function of a for fixed cj>. The parameter 
values are fixed as in Fig. [3] 



FIG. 7. (color online) Josephson currents in the p y -F-d yz 
junction, (a) Charge and (b) z-component spin currents as 
a function of <j> for fixed a. (c) Charge and (b) ^-component 
spin currents as a function of a for fixed <j>. The parameter 
values are fixed as in Fig. [3] 



Note that the spin current has opposite sign compared 
to the p z -wave junction as expected. 

In Fig. ^ we present the currents for the (p z + ip y )- 
wave TS state. As predicted in Sec. [TTJ the results for 
this junction are very similar to those for the p z -wave TS, 
as the lowest-order Josephson coupling proceeds through 
the p 2 -component of the chiral p-wave gap. The results 
are therefore summarized by Eqs. (53) and (54). 



Pz + 'Py " F - d v 



C. d yz -wave singlet superconductor 

The parity requirement for the superconducting gaps 
leads us to expect qualitatively different behaviour upon 
replacing the s-wave superconductor by a d yz -wave su- 
perconductor due to the even and odd dependence on k y , 
respectively. To test this, we repeat the above analysis 
for a dyz-w&ve pairing symmetry in the singlet supercon- 
ductor. Starting with the p z -w&ve TS, in Fig. © wc find 
that the charge current is approximately given by 



I c = I c sin(20) + I' c sin(20) cos(2a) , 



(55) 



with I c and I' c of comparable magnitude. This is clearly 
consistent with the predicted absence of single-Cooper- 
pair tunneling processes when only one of the order pa- 
rameters is odd in k y . The spin current is 0- independent 
and has the approximate form I$, z = I s sm(2a), charac- 
teristic of the contribution due to spin-flip reflection. 

In Fig. [7] we show the charge and spin currents for the 
p y -wave TS pairing symmetry. Since both the singlet 
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FIG. 8. (color online) Josephson currents in the (p z + ip y )-F- 
d yz junction, (a) Charge and (b) z-component spin currents 
as a function of 4> for fixed a. (c) Charge and (b) ^-component 
spin currents as a function of a for fixed <j>. The parameter 
values are fixed as in Fig. [3] 



and triplet gaps are odd in k y , we expect that a charge 
current is realized at the lowest order of perturbation 
theory. Indeed, as can be seen in panels (a) and (c), the 
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numerical results for the charge current are well approx- 
imated by Eq. (53), and the charge current is an order 
of magnitude larger than for the p z -wave junction. As 
expected, there is also a ^-dependent contribution to the 



spin current, which is again of the form Eq. (54) 



To conclude our survey, in Fig. [8] we plot the charge 
and spin currents for the chiral p-wave TS state. In 
this junction the d yz symmetry couples to ip y compo- 
nent of the (p z + ip y )-w&ve symmetry, and so the factor 
of i gives an additional ir/2 phase shift. Therefore, the 
lowest-order currents in this junction should be approx- 
imately obtained from the pure p y -wave case with the 
replacement <fr —> cj> + tt/2. Indeed the numerical results 
in Fig. [8] are consistent with the relations 



I c = I c sm((p) cos(a) , 



Is.z = Is sin(2a) + I s cos(y>) sin(a) 



(56) 
(57) 



IV. ROLE OF ANDREEV BOUND STATES 

The perturbation theory of Sec.[n]gives a good descrip- 
tion of the transport when multiple-Cooper-pair tunnel- 
ing processes make an insignificant contribution to the 
current, i.e. when higher-order terms in the perturba- 
tion expansion can be neglected. Although this is always 
the case sufficiently close to the transition temperature 
of the superconductors, resonant tunneling through An- 
dreev bound states can cause large deviations from per- 
turbation theory predictions at zero temperature. This 
is often particularly pronounced in junctions where the 
gap of the superconductors is odd in t he mo mentum com- 
ponent perpendicular to the interface! 28 ! ' It is therefore 
likely that in some of the junctions studied above, e.g. the 
junction between a p z -wave TS and an s-wave SS, there 
will be strong contributions to the current from higher 
harmonics in <j) and a at low temperatures. In order to es- 
timate the importance of this effect, here we consider an 
analytically-tractable model of the TFS junction where 
the current is due entirely to tunneling through Andreev 
bound states. 



This agrees with the predicted currents in Table [T] 

In contrast to the charge current, the spin current 
hardly changes upon replacing the s-wave SS by the d xy - 
wave SS. This reflects the dominance of reflection pro- 
cesses, which are insensitive to the superconductor on 
the other side of the junction. 



A. One-dimensional model junction 

We study a one-dimensional continuum model of the 
TFS junction with p z - and s-wave orbital symmetries for 
the TS and SS, respectively. Our analysis here closely 
follows that of Refs. Ql and El for a TFT junction. 
The energies of the Andreev bound states are obtained 
by means of solving the Bogoliubov-de Gennes equation 
H^{z) = E&(z), with Hamiltonian 



-e ia U M S(z) -&& + U S(z) 

-iQ{-z)A t d z -e(z)A s e^ 
\G(z)A s e^ »e(-2)A t £ 



For simplicity we assume that the effective mass m and 
chemical potential /z are the same on either side of the 
junction, and hence the Fermi wavevectors kF,v in each 
superconductor are also the same, i.e. kp v — kp. The TS 
and SS are described by the pairing potentials A t and A s , 
respectively. The barrier is modeled as a (^-function with 
charge scattering potential Uq and magnetic scattering 
potential Um = Um (cos ae x + sinae^). 

The Andreev bound states have energy lying within 
the bulk gap of the two superconductors, i.e. \E\ < 
mm{kp\A t \, |A S |}, and are hence exponentially local- 
ized at the interface. Within the Andreev approxima- 
tion, where the superconducting gap is assumed negligi- 
ble compared to the Fermi energy, an appropriate ansatz 



-iQ{-z)A t f z Q(z)A B e-** \ 

,i -Q(z)A s e-^ i0(-z)A t & 

e- ia U M S(z) - U S(z) + M / 



for the wavefunction of these states is 

¥ v (z) = e - VK " z (»f^ + e +ikFZ + *„ f _e- i * F *) . (59) 

Following the notation of Sec. [ll| v = s (t) denotes the 
singlet (triplet) side, and as a factor v = 1 (— 1). k„ is 
the inverse decay length in the v superconductor. The 
spinors ^ Vy ± are defined 

*t,± = (ot,±, b t , ± , ±e ±l T'a t , ± , T e±^6 t ,± ) T , (60a) 

*.,±=(o a ,±, b a ,±, -e«<*FT>&,,±, eW^'Vi ) T ' 

(60b) 



11 



where the subscript {+, — } indicates the direction of 
propagation, and the phases j v are given by 



cos 74 



cos 7 S 



E 
k F \A t \ 
E 



sin7( 



I A., 



Bin 7« = 



h 2 K t 

m\A t \ 
h 2 K s kp 



|A. S 



(61a) 
(61b) 



The a Vt ± and b v ^± appearing in Eq. (60) are constants to 
be determined by the boundary conditions obeyed by the 
wavefunction at the interface. In addition to continuity 
of the wavefunction across the junction, 

¥ t (z = 0-)=¥ a (* = + ), (62) 

we require that the derivative obeys 

0,*,(*)U=o+ -0**tC*)U=o- 

\ 

2k r '•' " ° °, 9,(z = + ), 

(63) 

in order to conserve probability. Here we use the dimen- 
sionless parameters 
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(64) 



to characterize the strength of barrier potentials. The 
boundary conditions give eight equations for the coeffi- 
cients a v .± and b u ,±, and have nontrivial solution when 
the energy of the bound state satisfies the equation 

= 2 + 8g 2 + 8g 2 cos(2a) - 2 cos(2</>) - W cos 2 {j t ) 

- [4 + 8g 2 + 8g 2 cos(2a)] cos 2 (7 s ) 
+ [4 + W + 16<? 2 ] cos 2 ( 7t ) cos 2 ( 7s ) 
-2 [2g 2 + 2Z 2 + l] sin(27 t ) sin(2 7s ) 
+16g cos(a) sin(</>) cos(7 t ) cos(7 s ) 

— 16<?cos(a) sin(</>) sin(7 t ) sin(7 s ) , 



where 



W = 4 (2Z 2 + l) + 16 (g 2 - Z 2 



(65) 



(66) 



B. Analytical solution at fei?At = A s 

The parameters that characterize triplet and singlet 
superconductors are, of course, independent from each 
other, as the two superconductors cannot be made from 
the same material. Therefore there are numerous choices 



to realize the parameters in Eq. (64), according to the 
materials that the junction is made. Nevertheless, in 
the limit when the pairing potentials of the singlet and 
triplet superconductor are the same, i.e., kpA t = A s , 
it is possible to express the Josephson charge and spin 



currents entirely in terms of the Andreev bound state 
energies. Although this is a highly-idealized situation, it 
clearly reveals the influence of resonant tunneling on the 
currents, which we expect to remain qualitatively valid 
for a more realistic model of the junction. 

When kpA t — A s , the bound state energy parameters 
are equal 7l = Jr, and Eq. (65) is drastically simplified 



1 

d 2 
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B = — cos(a) — - sin((/>) . 
2 x ' A 



(67) 

(68a) 
(68b) 
(68c) 



The positive Andreev bound state energies are hence 
found to be 

^ = VD\VDA + B± VDA - B\ . (69) 
\A S \ 

The Josephson charge and z-spin currents are definecP^' 

e ^ dEi Ei 

!c = - * > . «T tanh ^77^ > ( 70a ) 

Innb . (70b) 



1 



Is . Z , 



\ " 



BE, 



4 ^ da 

l—a,b 
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Inserting Eq. (|69|) into Eq. (|70|) we hence obtain 

Ir 
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DA + oB 



e|A s | x ^ 

8h ^-f 

x [2Dg cos(a) cos(</>) + a cos(0)] 
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(71a) 



a=± 



D 

DA + aB 



g 2 sin(2a) — Dg sin(a) sin(<^>) + ag sin(a) 



tanh ■ 



En 



2k B T 



crsign(_B) tanh 



E h 



2k B T 



(71b) 



We plot the Josephson currents at zero and finite temper- 
ature in the upper and lower panels of Fig. ([9]), respec- 
tively. To obtain the finite temperature results we assume 
that both gaps display BCS weak-coupling temperature- 
dependence with critical temperature T c . We find that 
the barrier potential Z mainly affects the amplitude but 
not the form of the current-phase relation, so hereafter 
we only present results for Z = 0. 
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FIG. 9. (a) Josephson charge and (b) spin current for the 
case g = 1, Z = 0, at zero temperature T = 0. (c) Josephson 
charge and (d) spin current for the same parameter set at 
finite temperature T/T c — 0.5. 




FIG. 10. (color online) (a) Josephson charge current Ic in 
the quasiclassical calculation at g = 1, Z = 0, T = 0. (b) Ex- 
pected behaviour of I c in the perturbative regime as given by 
Fi(a,(f>) = cos(a) cos(0) [Eq. \22)]. (c) Josephson spin cur- 
rent Js jZ in the quasiclassical calculation a,t g — 1, Z — 0, T = 
0. (d) Expected behaviour of Js.z in the perturbati ve r egime 
as given by F2(a,(j>) = sin(a) sin(</!>) — sin(2a) [Eq. ( 24 1] . 



The zero temperature results for the charge and spin 
currents at a = are in good agreement with the per- 
turbation theory predictions. Rotating the magnetic 
moment towards the y-z plane, however, we observe 
the current-phase relations display sharp jumps at the 
zero-energy crossings of the Andreev bound states. The 
Fourier decomposition with respect to the phase differ- 
ence hence contains a large contribution from higher har- 
monics, and thus multiple Cooper pair tunneling pro- 
cesses are important in the zero-temperature limit. In- 
deed, we see in Fig. ([9]) (a) that the maximum current at 
a = 0.57T is comparable to that at a = 0, whereas in 
the tunneling regime we expect that the latter should be 
much larger than the former. In contrast, the currents 
at half the critical temperature [panels (c) and (d)] are 
in much better agreement with the perturbation theory 
predictions and the lattice model calculations. In partic- 
ular, the amplitude of the charge current-phase relation 
at a = 0.57T is now much smaller than that at a — 0. 

To visualize the deviation from the perturbative re- 



sult at zero temperature, in Fig. ( 10 ) we plot side-by- 



side the exact and perturbative charge and spin currents 
as a function of angle a and phase (f>. Comparison of 
panels (a) and (b) clearly shows that the perturbation 
theory captures the broad features of the charge current. 
The greatest deviations occur near to a = ±0.57r, where 
we find jump discontinuities in the current due to zero- 
energy crossings of the Andreev bound states which are 
present for ||a| — 0.57r| < arcsin(l/2g). For g ^> 1, the 
zero-temperature current is well approximated by 



eA 1 
2h g* 



;ign(cos(a)) cos(4>) . 



(72) 



Note that the amplitude of the current goes as g 2 , in- 



stead of g J as predicted in Eq. ( 22 ): enhancement of the 



low-temperature current above the perturbation theory 
predictions is a well-known consequence of the presence 
of zero-energy Andreev statesP^We now turn to the spin 
current: comparing Fig. (10)(c) and (d), we again find 



that the perturbation theory reproduces the main out- 
lines of the exact Bogoliubov-de Gennes calculation. In 
particular, the functional form of spin current is best 
fitted by F 2 (a, </>) = — sin(a) sm(4>) — sin(2a), i.e. the 
contribution from sin(2a) term in Eq. ( |24[ ) is comparable 
to the sin(a) sin(</()) term. In the large-g limit, however, 
the phase-dependent part of the spin current becomes 
negligible, and we find 



Is.z = 



Aq 1 

" 4 g 



sign(cos(a)) sin(a) 



(73) 



Although the current is again enhanced beyond the per- 
turbation theory predictions, the two approaches agree 
that reflection processes dominate in the limit of a strong 
magnetic barrier. 

The critical charge current, defined as the maximum 
current with respect to <f>, is a readily-accessible experi- 
mental quantity. We plot this alongside the critical spin 



current in Fig. ( 11 ) as a function of the temperature and 



the angle a. As can be seen in panel (a), the temperature- 
dependence of the critical charge current qualitatively 
changes with the orientation of the magnetization: For 
a = it grows linearly with decreasing temperature im- 
mediately below T c , and saturates at T ss 0.2T C ; In con- 
trast, when a = 0.57T, it grows superlinearly with decreas- 
ing T, but remains much smaller than the a = current 
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FIG. 11. (a) Critical Josephson charge and (b) spin current as 
a function of temperature, at several different values of a. (c) 
Critical Josephson charge and (d) spin current as a function 
of a, at several different temperatures. In all panels we take 
the magnetic potential g = 1. 



until T 0.2T C , below which it displays a rapid increase. 
This "low-temperature anomaly" in the a = 0.5-7T crit- 
ical charge current reflects the importance of resonant 
tunneling through the zero-energy Andreev bound states 
far below T C P^ Similarly, in panel (c) we observe that 
near T c the critical current as a function of a follows the 
perturbation predictions of max{|/ c |} oc | cos(a)|; as the 
temperature is decreased towards T = 0, however, the 
critical current at a = ±0.57r increases due to the reso- 
nant tunneling, and there is hence overall relatively weak 
dependence of the critical current on the magnetization 
orientation. The critical spin current [Fig. (ll)(b), (c)] 
also shows strong temperature- and a-dependence. 



CONCLUSIONS 



In this paper we have studied the unconventional 
Josephson charge and spin currents in a TFS junction. 
Using complementary theoretical methods, we have es- 
tablished that single-Cooper-pair tunneling currents are 
possible when the magnetization of the FM has a com- 
ponent parallel to the d vector of the TS, and when the 
orbital pairing states of the superconductors do not have 
the same parity with respect to the interface momentum. 
We hence see that spin and orbital degrees of freedom 
both play a critical role in this junction; this is also the 
case when a lowest-order Josephson effect between a TS 
and an SS is mediated by spin-orbit coupling at the bar- 
r j er J30H33] ^ a microscopic level, the perturbation theory 



analysis reveals that the spin-dependent phase shifts of 
the tunneling Cooper pairs are responsible for the charge 
and spin Josephson currents, due to the interference of 
the spin-f and spin- 4- particle currents in the TS. Sur- 
prisingly, the interference of the particle currents is also 
responsible for a phase-dependent spin current in the TS, 
even though spin currents are forbidden in the SS. Simi- 
lar interference effects occur in the TFT junction ,^ an d 
junctions between TSs with misaligned d-vectors! 14 * 15 ! 

Our analysis has confirmed previous observations of a 
highly unusual charge current in the TFS junctionPSETJ 
Not only is it linearly proportional to the magnetization 
of the FM as oc M • d, it also has unconventional cos<?!> 
dependence on the phase difference (if time-reversal- 
symmetry is not broken in one of the superconductors) . 
This implies a contribution to the junction free energy 



M • dA,A t sin< 



(74) 



In this expression we can regard dA s A t sin0 as an in- 
trinsic interface spin which appears at junctions between 
a TS and an SS, even when the barrier is nonmagnetic. 
If the barrier is ferromagnetic, the coupling of its mag- 
netic moment to this intrinsic spin therefore generates 
the lowest-order Josephson effect. Remarkably, such an 
intrinsic interface spin is i ndeed kno wn to exist in non- 
magnetic TS-SS junctions.^ 29 ' 49 ' 50 ' For an s-wave SS, 
this spin only appears for exactly the same p-w&ve TS 
orbital configurations which would allow a lowest-order 
Josephson current in the corresponding TFS junction. 

The unusual form of the Josephson current in the TFS 
offers strong tests for a triplet state. For instance, the 
observation of the linear dependence of the current on 
M would be clear evidence of triplet pairing. On the 
other hand, a domain structure in the ferromagnet could 
significantly reduce the Josephson current, as the cur- 
rents across domains with opposite magnetization would 
have opposite sign. This can be turned to our advantage, 
however, as a magnetic flux line trapped at the bound- 
ary between two such domains would be quantized in 
half-integer multiples of <E>o- This is a key signature of 
the lowest-order Josephson coupling in the TFS junction, 
and could be directly imaged with SQUID microscopy, 
or deduced from the Fraunhofer pattern. More specula- 



tively, the coupling Eq. ( 74 1 could spontaneously induce a 



magnetization in a barrier sufficiently close to a magnetic 
instability, if the free energy gain due to the Joseph son 
coupling can offset the cost of magnetic energyP^H 

In our study we have neglected the likely variation of 
the superconducting order parameter close to the junc- 
tion interface. Since our results depend only on the bulk 
properties of the superconductors, however, we do not 
expect qualitative modification of our results. A more 
serious limitation of our calculation is that we have not 
accounted for the torque exerted by the spin current on 
the barrier's magnetic moment. Regarding the d-vector 
as fixed, we anticipate that the polarization of the spin 
current oc d x M would cause a precession of the mag- 
netization about the d-vector, with eventual decay into 



14 



the stable configuration.^ If the d- vector is only weakly 
pinned, on the other hand, there may be a significant 
reconstruction of the TS pairing state close to the inter- 
face. Although this is a very interesting problem, it is 
beyond the scope of the current paper. 
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